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ABSTRACT 

Approximately half the baryons in the local Universe are thought to reside in the warm- 
hot intergalactic medium (WHIM), i.e. diffuse gas with temperatures in the range 
10^ K < T < 10^ K. Emission lines from metals in the UV band are excellent tracers of 
the cooler fraction of this gas, with T < 10^ K. We present predictions for the surface 
brightness of a sample of UV lines that could potentially be observed by the next 
generation of UV telescopes at z < 1. We use a subset of simulations from the OWLS 
project to create emission maps and to investigate the effect of varying the physical 
prescriptions for star formation, supernova and AGN feedback, chemodynamics and 
radiative cooling. Most models agree with each other to within a factor of a few, 
indicating that the predictions are robust. Of the lines we consider, C III (977 A) is 
the strongest line, but it typically traces gas colder than 10^ K. The same is true for 
Si IV (1393,1403 A). The second strongest line, CIV (1548,1551 A), traces circum- 
galactic gas with T - lO^K. OVI (1032,1038 A) and NeVIII (770,780 A) probe the 
warmer (T ~ lO'^'^ K and T ~ 10^ K, respectively) and more diffuse gas that may be a 
better tracer of the large scale structure. N V (1239,1243 A) emission is intermediate 
between C IV and O VI. The intensity of all emission lines increases strongly with gas 
density and metallicity, and for the bright emission it is tightly correlated with the 
temperature for which the line emissivity is highest. In particular, the C III, C IV, 
Si IV and O VI emission that is sufficiently bright to be potentially detectable in the 
near future (surface brightness > 10'^ photon cm~^ sr~^), comes from relatively 
dense (p > lO^pmcan) and metal rich {Z > O.IZq) gas. As such, emission lines are 
highly biased tracers of the missing baryons and are not an optimal tool to close the 
baryon budget. However, they do provide a powerful means to detect the gas cooling 
onto or flowing out of galaxies and groups. 

Key words: method: numerical - intergalactic medium - diffuse radiation - radiation 
mechanisms: thermal - cosmology: theory - galaxies: formation 



1 INTRODUCTION 

In the high redshift Universe most baryons are believed 
to reside in gas that is traced by the Lyg forest which is 
seen in the spectra of quasars (e.g. iProchaska &: TumlinsonI 
I2OO8I for a review). In the low redshift Universe, how- 
ever, the LyQ forest may trace less than half the baryons. 
Numerical simulations suggest that this is due to the 
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progressive heating of the diffuse intergalactic medium 
(IGM) by gravitational shocks produced when structures 
form (e.g. fsunvaev fc ZerdovichI Il972l: iNath fc siilj I2OOI 



Furlanett o fc Loebl 12004 iRasera fc Tevssieij l2006l : iMeiksinl 
2009 for a review). 

Approximately half of the intergalactic gas is predicted 
to have temperat ures in the range 10^ K < T < 10^ 
K at 2 » (e.g. ICen fc Ostrikcr 19^; iDave et all I2OOII : 



iBertone. Schave fc Doiad[2008i for a review). As yet, a large 
fraction of this shock-heated gas, or warm-hot intergalactic 
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medium (WHIM, hereafter), has not been unambiguously 
detected. 

Metal line transitions in the ultraviolet (UV) band pro- 
vide a promising route to detecting the cooler fraction of the 
WHIM with 10^ K< r < 10" K. The OVI doublet is ideal 
for tracing this whole temperature range, and a handful of 
other transitions may help to identify gas in either slightly 
cooler or warmer temperature regimes. In particular, if col- 
lisional ionisation dominates, as we will show to be the case 
for gas that is dense enough to be detectable in emission, 
the CIV and N V doublets trace gas with T ~ 10^ K, while 
the NeVIII doublet is produce d by gas with T ~ 10*^ K 
(e.g. [Sutherland fc Dopita|[T99l l . The hottest fraction of the 
WHIM, with 10' ^ K< T < W'^ K , is better tr aced by X-ray 
trans itions (e.g. lBregmanll2007l for a review; iBertone et al.l 
l2010t ). 

Cooler WHIM gas (T < lO" K) may have been suc- 
cessfully observed in absorption in FUSE and STIS spectra 
in the UV band, although the interpretation of these 
observations remains controversial. Narrow Lya, C IV and 
O VI absorption lines in QSO spectra hav e been detected in 
large number s in hi g h resolution spectr a ( Tripp et al. 2000l: 

33; 



Savage et al. ' '2002"; toanforth & Shull' '200- 
20081 : iDanforth & ShuU 200^; ,Cooksev et al. 



Tripp et all 



E p et al.l 
l2008l: 



Richter et al] I2OO8I: Pfhom fc ChenI l2008al: iThom fc ChenI 



2008bl : IWakker fc Savagel l2009l : ICooksev et al.l bOld V In 
tergalactic NeVIII, w hich, so far, has bee n observed only 
in the Galactic halo fea vage et all l2005ll , and NV will 
also potentially be detected by the Cosmic Origin Spec- 
trograph ( COS. iFroning fc Greenll20oi : iPaerels et al]|2008l : 
iRichter et al.l |2008| ) on board the Hubble Space Telescope 
(HST). 

While narrow absorption lines are believed to trace 
relatively cool and predominantly photo-ionised gas 
iHowk et al.ll2009l ). broad lines have been interpreted as the 
signature of mostly coUisionally ionised gas with tempera- 
tures in excess of 1 0° K iRicht cr ct al. 2004; Richter ct al] 
20061: iLehner et al .' 2007; Wakkcr fc Savag^ |2009| . but see 
Oppenheimer fc Dave ,2009. ). The detection of broad Lya 
absorption lines with velocity widths larger than b > 40 
km s~^ in FUSE and STIS sp ectra (Richter et al.l |2004|: 
Sembach et"ai] l2004l: iLehner et al. 2007: Danforth fc Shulj 



20ld) further improves our understanding of the WHIM and 



opens the way to accounting for an even larger fraction of 
the missing baryons and to constrain the sources of the IGM 
metal enrichment. 

While UV absorption lines are a powerful method by 
which the physical properties of the WHIM can be probed, 
they only provide 1-dimensional information along the line 
of sight. Multiple lines of sight from close quasar pairs 
are necessary to reconstruct the 3-dimensional gas distri- 
bution, a challenging task that may be within the reach 
of COS. Conversely, emission lines are able to efficiently 
probe the full 3-dimensional matter distribution by map- 
ping the 2-dimensional distribution of the gas on the sky 
and by adding information about the thi rd dimension by 
means of the spectral reds hift of the lines jFurlanetto et all 
I2OO4I : ISembach et alll2009l ) . Because emissivity scales as the 
gas density squared, emission lines are biased towards prob- 
ing higher density regions than absorption lines. Therefore, 
while emission may provide an ideal probe of the proper- 
ties of the high density WHIM, absorption lines may re- 



Table 1. List of emission lines. The first column shows the ion 
giving rise to the line, the second and the third columns the wave- 
lengths of the two components of the doublet. The strongest com- 
ponent is listed in column 2 and its energy in column 4. C III is 
the only singlet in the sample. 



Ion 


Al (A) 


A2 (A) 


El (eV) 


cm 


977.03 




12.690 


CIV 


1548.187 


1550.774 


8.008 


Si IV 


1393.755 


1402.770 


8.896 


NV 


1238.821 


1242.804 


10.008 


OVI 


1031.912 


1037.613 


12.015 


NeVIII 


770.409 


780.324 


16.094 



main the most powerful tool when in vestigating mildly dense 
and underdense regions (see e.g. iFurlanetto et al.l |2004| : 
iBertone. Schave fc Dolaell2008l : iBertone et al.ll2010l l. 

In this work we employ a subset of hydrodynamical, 
cosmological simulations from the Ove rwhelmingly Large 
Simulations (OWLS, hereafter) project l|Schave et al.|[2O10h 
to predict the intensity of a sample of UV emission lines that 
could potentially be detected by current and future instru- 
ments, such as the Fa int Inte rgalactic Redshifted Emission 
Balloon ( f/j^^^^ Z/Z/. iTuttle et al. 2008.) and the Advanced 
Technology Large A perture Space Telescop jj] (ATLAST 
Postman et al.|[2008l ). In a companion paper (|Bertone et al.l 



2010l . Paper I in the following) we have investigated the soft 



X-ray emission from the WHIM using the same methodol- 
ogy, and we refer to that work for more details about the 
numerical aspects of our calculations. 

The OWLS runs are ideal for this study because they 
couple a large simulated dynamical range at high resolution 
with a large set of models in which the implementation of 
several physical processes has been varied. This feature in 
particular allows us to investigate the dependence of the 
predicted UV emission on a number of uncertain physical 
prescriptions, in addition to considering standard tests for 
mass resolution and the size of the simulation box. In this 
work we use simulations of 100 comoving Mpc on a side 
and focus on the strongest line of the UV doublets C IV, 
Si IV, NV, O VI and NeVIII and the singlet CIII, listed in 
Table [1] The intensity of the weaker line in the doublets can 
be recovered by multiplying the intensities derived in this 
study by a factor of 0.5, which corresponds to the difference 
in the line emissivities in the doublet. 

This paper is organised as follows. We briefly describe 
the simulations and the method by which we calculate the 
gas metal line emission in Section 12.11 Section [3] contains 
our results for z = 0.25, while in Section |4] we test the effect 
of varying the angular resolution used to build the maps. 
In Section Owe consider the redshift dependence of the flux 
intensity. Section |6] investigates what kind of gas produces 
most of the emission and Section [7] examines the depen- 
dence of the result s on the physical model . We compare our 
results to those of IFurlanetto et al] l|2004h in Section [8] and 
we discuss the prospects for observing the UV emission with 
future space-borne telescopes in Section |9l Finally, we sum- 
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marise our conclusions in Section [TO] and describe numerical 
convergence tests in the Appendix. 

Throughout thi s paper we will a ssume a WMAP3 
ACDM cosmology jSpergel et al.l |2007| ) with parameters 

= 0.238, f^b = 0.0418, Qa = 0.762, n = 0.951, and 
(Tg — 0.74. The Hubble constant is parametrised as Hq = 100 
h-^ km s"^ Mpc"\ with h = 0.73. 



2 NUMERICAL METHODS 



Table 2. Adopted solar abundanc es, f rom lAI Iende P rieto et ahl 
ll200ll' l. lAllende Prieto et al.l bOO'A and lHolweged 1 I2OOII ). 

Element ni/nn Element rii/nn 



In Section [2.11 we briefly describe the subset of simulations 
from the OWLS project employed in this study, while in 
Sections 12.21 and 12.31 we review the methods by which we 
calculate emission line intensities, and generate synthetic 
surface brightness maps, respectively. 



2.1 Numerical simulations 

In this work we employ a subset of the cosmological, hy- 
drod ynamical simulations that comprise the OWLS project 
IScha vc ct al. 2010). A detailed description of the numer- 
ical methods emplo yed in the OWLS project is given in 
ISchave et al.l (|2010l ). Here we give a brief overview of the 
subgrid physical modules used in the reference model (REF, 
hereafter) that is used to derive most of the results in this 
paper. In Section[7]we discuss a number of other models that 
incorporate variations of the physical modules, implemented 
one at a time. 

The simulations were run using a significantly ex- 
tend ed version of the parallel PMTree-SPH code GADGET 
III (|Springel|[2005l ). Here we use runs with box sizes 100 

comoving Mpc containing 512^ particles of both gas 
and dark matter. The baryonic particle mass is 8.66 x 10^ 
h-^ Mq. 

Gas cooling i s implemented following the prescription 
of|w lersma et al.l (,2009a fl. Net radiative cooling rates are 
computed and tabulated element by element for all 11 
elements (H, He, C, N, O, Ne, Mg, Si, S, Ca and Fe) 
tracked by the simulation in th e presence of the cosmic 
microwave background and the iHaardt fc Madaul (|200ll ) 
model for the UV/X-ray background radiation from quasars 
and galaxies. The tables, created with the publicly avail- 
able photo-ion ization package CLOUDY (last described by 
iFerland et all 1998) assume the gas to be optically thin and 
in (photo-)ionization equilibrium, are interpolated over den- 
sity, temperature and redshift. 

Star forma tion is implemented following the sub-grid 
prescription of ISchave fc Dalla Vecchial (|2008l ). Gas with 
densities exceeding the threshold for the onset of therm o- 
gravitational instability (nu = 0.1 cm~^: [Schavd |2004| ) is 
considered to be multiphase and follows an effective equation 
of state of the form P oc p'^"'^ , normalised to P/k = 1.08x 10^ 
K cm~^ at the density threshold. Multiphase gas forms stars 
at a pressure-dependent rate that reprod uces the observed 
Kennicutt-Schmidt law (jKennicuttll 19981 ). The prescription 
by which the simulations track the subsequent chemical en - 
richment of the gas is described in .Wiersma et al.l (|2009bl ). 



^ Using their equation (3) rather than (4) and CLOUDY version 
05.07 rather than 07.02, as used in that paper. 



H 

He 

C 

N 

O 

Ne 



1 

0.1 

2.46x10-4 
8.51x10-5 
4.90x10-"' 
1.00x10-4 



Mg 

Si 

S 

Ca 
Fe 



3.47x10-5 
3.47x10-5 
1.86x10-5 
2.29x10-'' 
2.82x10-4 



In brief, we follow the timed release of 11 different ele- 
ments from massive stars (Type II supernovae and stellar 
winds) and intermediate mass stars (Type la supernovae 
and a symptotic giant branch stars), assuming a Chabrier 
l|Chabrier.,2003. ) IMF in the mass range 0.1-100 Mq. 

Energy injection from supernovae is included 
as kinetic feedback f ollowi ng the prescription of 
iDalla Vecchia fc Schavd (|2008h . Two parameters are 
used to describe the energy injection: the wind velocity, 
Hw, and the mass loss rate, 77 — rh^/rh,, which describes 
the amount of gas kicked into the wind, a function 

of the rate of star- formation, m*. The REF model uses 
r) = 2 and «„ = 600 km s~^, which correspond to 40% of 
the supernova energy being available in the form of kinetic 
energy. The REF model does not include a prescription for 
AGN feedback, but we will consider a model with AGN in 
Section [T] 



2.2 The emissivity tables 

The intensity of the emission lines is calculated following the 
procedures of Paper I. We generated tables of the line emis- 
sivity as a function of hydrogen number density, temperature 
and redshift for about 2000 of the strongest lines of the 11 
elements tracked by OWLS. The line e missivities are calcu- 
lated using CLOUDY (version c07.02.02. [Perland et allll998l . 
released in July 2008) for an optically thin gas in ionisation 
equilibrium in the presence of the cosmic microwave back- 
groun d and a uniform, evolvin g meta-galactic UV radiation 
field (|Haardt fc Madaul I2OO1I ) from galaxies and quasars. 
The same assumptions were used in the calcu lation of the 
radiative cooling rates (| Wiersma et al.l|2009al ). The tables 
sample the temperature range 10^ K < T < 10*'^ K in bins 
of ALogj^yT = 0.05 and the hydrogen number density range 
10"* cm-^ < wh < 10 cm"^ in bins of ALogignn = 0.2. 

The tables were computed for solar metallicity, with 
the adopted solar abundance Zq = 0.0127. This cor- 
responds to the value obtained using the default abun- 
dance set of CLOUDY, listed in Table [21 which is a 
combi n ation of the abund a nces of lAllende Prieto et al.l 
(|200ll ). lAllende Prieto et al.l (|2002l ) and iHolwegerl l|200ll ). 
Our adopted abu ndance set ma y differ significantly from 
that estimated bv lLoddersI (|2003l ). In CLOUDY, the oxygen 
abundance in particular is about 20 per cent lower than es- 
timated bv lLodderd l|2003D . The differences in the assumed 
element abundances should be kept in mind when comparing 
with predictions from different codes. 

Fig. H] shows the normalised emissivities of the UV lines 
listed in Table [T] as a function of temperature, for constant 
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Figure 1. The normalised emissivity, e/n^, in units of erg cm^ s~^, of a selection of UV emission lines as a function of temperature 
assuming solar abundances at z = 0.27. The left, middle, and right panels show results for constant hydrogen number densities rtfj = 1, 
10-3 and 10"*^ cm respectively. If the gas is coUisionally ionised, the normalised emissivity is independent of the gas density. Cool 
diffuse gas is mostly photo-ionised, while collisional ionisation equilibrium dominates for the highest temperatures and densities. 



hydrogen number densities nu 



1, lO-'' and lO"" cm" 



(left, middle and right panels, respectively). Fig. [T] demon- 
strates that collisional ionisation dominates at the highest 
temperatures and densities, while photo-ionisation becomes 
the main ionisation process only when both the density and 
the temperature are very low. This can be seen by compar- 
ing the three panels of Fig. [T] and noting that the normalised 
emissivity peaks at a well defined temperature and is inde- 
pendent of the density, if collisional ionisation dominates. 
As we will see in the following, most of the gas emission is 
produced in coUisionally ionised regions with gas densities 
nu > 10-3 cm-3. 

The assumption of ionisation equilibrium, a l so use d 
to calculate the cooling tables (jWiersma et al.l l2009al \ 
is justified for photo-ionised reg ions and dense g as in 



the centres of clusters (see iBertone. Schave fc Dolae |2008| 
for a review). However, non-equilibrium processes may 
become important in the WHI M and in the outer 
regions of groups and clusters |Hughes &: Singhl 19941: 
Yoshida et al. 2005: Yoshikawa & Sasaki"2006': 'Cen & FanJ 
2006 ; Gnat fc Sternber g 2007 ) . Gnat fc Sternberg (2009. ) , in 
particular, have recently shown that non-equilibrium pro- 
cesses affect the cooling of shock-heated gas with temper- 
atures 10* K < T < 10^ K, and that the amplitude of the 
deviations from ionisation equili brium incre ases with t he gas 
metallicity. On the other hand, lYoshikawa fc Sasa kl Hooi) 
found the effect of non-equilibrium ionisation to be small for 
WHIM emission. Finally, we note that because these studies 
ignored photo-ionisation, they may have overestimated the 
importance of non-equilibrium ionisation. 



2.3 The metal line emission 

The emission of the sample of UV lines listed in Table [T] 
is calculated following the procedures described in Paper I 
and we refer the interested reader to that work for the full 
details. 



We calculate emission only for non-star forming gas 
with TiH < 0.1 cm"3 because our simulations lack the res- 
olution and the physics to model the emission from higher 
density gas in the interstellar medi um, which is expected to 
be multiphase (see the discussion in lSchave fc Dalla Vecchial 
2008). 

The luminosity of emission line I for gas particle i, in 
units of erg , is calculated as a function of its gas mass 



density pi, hydrogen number density nH,i and ele- 



Li.i (z) = Eij^Q {z,Ti,n-a) 



X. 



y0 



(1) 



where £i,!,0 {z,Ti,nH) the emissivity of the line I, in units 
of erg cm-3 s^^, bi-linearly interpolated (in logarithmic 
space) from the emissivity tables as a function of the par- 
ticle temperature Logj^gTi and hydrogen number density 
Log]^gnH,i at the desired redshift. Xy.i is the mass frac- 
tion of element y and X^q the corresponding solar value. 
We use the "sm oothed" element abundances described in 
|w iersma et al.l (2009b), which alleviates the effect of the 
lack of metal mixing inherent to SPH simulations, although 
it does not solve the problem. The main effect of the use of 
smoothed abundances is to spread the emission over slightly 
greater gas masses, particularly in regions with lower metal- 
licities. 

The corresponding particle flux Fi^i, in units of photon 
s"^ cm"^, is given by 



-(1 + ^): 



(2) 



where _Dl the luminosity distance, hp the Planck constant, 
c the speed of light and Xi the rest-frame wavelength of the 
emission line I. 

To create maps of the line emission, the particle fluxes 
are projected on to a 2-dimensional grid using a flux- 
conserving smoothed particle hydrodynamics interpolation 
scheme. Finally, surface brightness (Sb, hereafter) maps are 
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Figure 5. The surface brightness PDFs of our sample of UV 
emission lines, also shown in Fig. \3\ and listed in Table [T] at 2; = 
0.25. The main plotting area shows only the high flux tail of the 
distributions, while the full PDFs are shown in the inset. The 
upper axis shows the flux per pixel in units of photon s~^ cm^'^. 
The pixel size is 15" (which corresponds to a physical size of 
42 kpc) and PDFs are calculated using five slices through the 
simulation box that are each 20 comoving Mpc thick (which 
corresponds to 1772 km s~^). In Appendix IA3I we show that for 
sufficiently large fluxes (> 10"'^ photon s~^ cm~^ sr~^ for O VI) 
the PDF is proportional to the slice thickness. The C III and 
C IV lines are the strongest emission lines in the sample, closely 
followed by OVI. NV and NeVIII fluxes are weaker than CIV 
by about an order of magnitude. 



computed by dividing the flux in each pixel, p, by the sohd 
angle Q it subtends: Ssp.i = Xii Fi,i/^- 

The number of pixels in the grid is determined by the 
angular resolution required to match the capabilities of cur- 
rent or future instruments. We investigate how fluxes change 
with grid angular resolution in Sectional A region with co- 
moving size 100 Mpc a.t z — 0.25 corresponds to a field 
of view of about 8 degrees for our chosen cosmology. An an- 
gular resolution of 15" on the sky, used throughout this work 
unless otherwise noted, requires a grid with 1926^ pixels. At 
z — 0.25 this angular resolution corresponds to a physical 
scale of 42 h"^ kpc, which is about 21 times larger than the 
maxi mum gravitat ional softening scale used in the simula- 
tion jSchave et al. 2010). We divide each simulated box of 
100 Mpc into 5 slices, each 20 Mpc thick, and cre- 
ate as many emission maps, which we combine to calculate 
the flux probability distribution function (PDF, hereafter) 
corresponding to this angular resolution and slice thickness. 
At z = 0.25, a slice thickness of 20 Mpc corresponds to 
an energy resolution AE/E ~ 6 x 10""^, which is roughly 
comparable to that of FIREBALL. 



Figure 6. As Fig. (5] but with the emission line flux expressed in 
erg cm~^ sr~^, instead of photon cm~^ sr~^. The relative 
strength of the lines varies slightly when converted to these units 
and is more easily comparable to the flux limits proposed for 
future missions. 

3 MAPS OF THE UV EMISSION 

In this Section we present emission maps for the strongest 
and most interesting UV emission lines we have identified 
as tracers of the cooler WHIM component with T < 10*^ K. 

Fig. [2] gives an impression of the morphology of the UV 
emission from large-scale structure. The left panel depicts 
the full simulation box and the two companion panels show 
two consecutive zooms into the core of a galaxy group. The 
displayed regions decrease in size by a factor of ten, ranging 
from 100 h-^ Mpc to 1 h'^ Mpc. 

Fig-Oshows surface brightness maps at z = 0.25 for the 
six UV lines considered in this work. The surface brightness 
PDFs are presented in Fig. [5] in units of photon s~^ cm~^ 
sr~^, and in Fig.[S]in units of erg s~^ cm"'^ sr~^. For compar- 
ison, the density, temperature, and metallicity in this region 
are shown in Fig. |4l which is the same as Fig. 13 of Paper 
I. This is the same region of space shown in all the emission 
maps in this paper, with the exception of Fig. [2l 

Fig. [3] reveals that none of the UV lines gives rise to 
much emission from the largest haloes (cf. Fig. [J) as in these 
haloes gas is shock heated to T >> 10^ K, far above the tem- 
peratures at which the ions responsible for UV emission are 
abundant. The OVI and NeVIII fluxes show similar spa- 
tial distributions with most of the emission associated with 
filaments. The CIII, CIV and Si IV emission maps demon- 
strate a much more clumpy morphology and their emission 
is associated chiefly with circum-galactic gas. This is a di- 
rect consequence of the fact that the peak of the emissivity 
shifts to higher temperature with increasing atomic number 
and ionization state, as can be seen from Fig. \T\ 

By comparing these results to those of Paper I, we 
find that the UV emission has a different spatial distribu- 
tion from the soft X-ray emission. While the X-ray emission 
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[Ht] 

6 1032A z = 0.25 




-4-2 2 

Log,o Sg (photon cm"' s"' sr"') 



Figure 2. Zoom into an OVI 1032 A emission map of a high-density region at z = 0.25. The object at the centre of the displayed 
region is a relatively small group. The slice thickness is 20 comoving Mpc (which corresponds to 1772 km s~^ for all panels). From 
left-to-right, the pixel sizes are 100", 10" and 1" and the angular sizes of the regions shown are 8 degrees, 48' and 4.8' (corresponding to 
100, 10, and 1 comoving Mpc), respectively. The blobs in the right panel are resolved substructures. 



mostly comes from the warmest gas in the largest haloes and 
filaments, UV lines trace smaller structures, such as galaxies 
and small groups. This difference can be partially explained 
by the different peak emissivity temperatures of the differ- 
ent ions. Here we find that OVI and NeVIII lines trace 
primarily the highest temperature, shocked gas (~ 10^ K) 
in filaments, in the surroundings of groups and in the in- 
tragroup medium itself, whereas GUI, CIV and Si IV lines 
trace gas at lower temperatures, have very inhomogeneous 
distributions and trace gas associated with galaxies. Si IV 
lines trace only very compact regions of space, that can be 
identified with gas in the haloes of galaxies. The spatial dis- 
tribution of C III emission is intermediate between those of 
Si IV and C IV. 

The PDFs in Figs. [5] and [6] demonstrate that C III is 
the strongest of the emission lines considered in this study, 
and the strongest fluxes are concentrated in very few pix- 
els (or equivalently, in very small regions). This is partly a 
consequence of the C III line emissivity (Fig. [T]) which, at 
high density, has a peak value at T ~ 10^ K that is larger 
than that of all other UV lines by 0.5 dex, and partly of the 
large amount of metals in the C'^^ phase, as we show in Sec- 
tion 16.11 C IV surface brightnesses are typically an order of 
magnitude lower than those of C III, in high density regions. 
N V emission is weaker than CIV by a factor of 10-100 and 
NeVIII by a factor of about 10. On average, OVI emission 
is fainter than C IV by a factor of about 2. The different 
spatial distributions of the ions that trace high tempera- 
ture gas (OVI, NeVIII, NV) and those that trace cooler 
gas (CIII, CIV, Si IV) mean that in spite of not being the 
strongest line in high density regions, O VI emission is able 
to exceed the emission from all other lines in low density 
regions, where carbon and silicon emission declines steeply. 
We therefore conclude that the best line in the UV band for 



probing the properties of the diffuse WHIM is O VI, whereas 
the strong lines such as C III and C IV that trace cooler gas 
offer the exciting opportunity of observationally probing the 
colder gas around, and accreting onto, galaxies. 

Fig. [7| shows the relative intensity of lines with respect 
to each other. Where contours lie above the dashed diago- 
nal line, the emission line on the i/-axis produces the high- 
est flux, when below, the line on the a;-axis is strongest. 
Line pairs with similar peak temperatures produce narrow 
contours that may lie on either side of the diagonal, while 
line pairs with different peak temperatures produce broad 
profiles. The large dispersion indicates that different lines 
are dominant in different density and temperature regimes. 
Fig. [7] illustrates the distribution of line ratios that can be 
expected when surveying the WHIM emission in a large field 
and may ultimately help to identify lines in different regions. 



4 ANGULAR RESOLUTION 

In this Section we test the effect of varying the angular res- 
olution used to make the emission maps. Fig. [8] compares 
the PDFs of the surface brightness in maps with resolu- 
tions i9 = 5" , 15" , 1', 4' and 10', which correspond to physical 
scales varying from 14 kpc to 1.661 Mpc. 

We find that most distributions, with the exception of 
the •& — 10' case, converge well at intermediate and low 
surface brightness. We find no convergence at the bright 
end, where the highest predicted surface brightness steadily 
increases with increasing resolution, at least on the scales 
considered here. This is a consequence of the fact that large 
pixel sizes prevent the detection of very strong emission com- 
ing from structures with an angular size smaller than that 
used to observe them. This said, diffuse gas on scales larger 
than galactic scales, such as the WHIM gas, can be observed 
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Figure 3. Maps of the emission for a sample of six different UV lines, listed in table[T] at z = 0.25. The figure shows a zoom of a region, 
with a comoving size of f4 Mpc (f .12 degrees on the sky), that includes the largest group in the simulation, with a total mass of a 
few times 10^* ^0- "^^^ density, temperature, and metallicity in this region are shown in Fig. |4] All panels assume the same angular 
resolution (15", which corresponds to a physical size of 42 kpc) and slice thickness (20 comoving Mpc, which corresponds to 
1772 km s-l). The maximum of the colour scale has been set to a fraction of the real maximum to enhance the emission of the weaker 
lines and of low-density regions. While strong CIV emission is concentrated around the haloes of galaxies, O VI and NeVIII lines trace 
more diffuse gas in a variety of structures such as filaments and haloes. Emission from UV lines is suppressed at the centres of large 
haloes, where the gas is hotter than 10'' K and X-ray emission dominates. 



with an angular resolution of about 10-15" at z ~ 0.25 with- 
out significant loss of information. Higher angular resolution, 
however, might be desirable when observing gas with more 
clumpy spatial distribution on smaller scales, such as out- 
flowing gas in galactic winds. For reference, in the OWLS 
cosmology a pixel size of 15" subtends a comoving (proper) 
size of about 52 (42) h'^ kpc at z = 0.25, 98 (66) ft"^ kpc 
at z = 0.5 and 175 (88) kpc mt z = 1. 



5 EVOLUTION 

In this Section we briefly discuss how the intensity of the 
emitted flux varies with cosmic time. 



In Fig. [9] we show the surface brightness PDF of O VI 
emission at z — 0.25, 0.5 and 1. The overall shape of the 
PDF does not vary substantially, but the highest predicted 
flux increases with decreasing redshift. This is at least partly 
due to the fact that at the constant angular resolution of 15" 
assumed for the maps, the physical size that corresponds to 
a pixel increases steadily with redshift, from 42 kpc at 
z = 0.25 to 88 kpc at z = 1. As shown in Section|4l the 
maximum flux decreases with increasing pixel size. However, 
the emissivity also depends on the density, temperature and 
metallicity of the emitting gas, all of which evolve with cos- 
mic time. 
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Figure 4. As Fig. [S] but siiowing maps of tiie gas density (left panel), temperature (middle panel) and metallicity (rigiit panels) for 
both star forming and non-star forming gas. 
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Figure 8. As Fig. [S] but for OVI emission as a function of the 
angular resolution of the emission maps. The pixel sizes are i9 = 
5", 15", 1', 4', and 10', corresponding to physical sizes of 14 
kpc, 42 h-^ kpc, 166 h'^ kpc, 665 h'^ kpc and 1.661 h'^ Mpc, 
respectively. With the exception of the i? = 10' and 4' cases, the 
results are converged at intermediate and low surface brightness 
(Sb < 10^ photon s'^ cm ^ sr ^). The maximum predicted 
surface brightness increases steadily with angular resolution, a 
consequence of the fact that very strong emission on scales smaller 
than the resolution is smoothed out. There is no clear convergence 
at low angular resolutions (i? > 1'). 



Figure 9. As Fig. [5] but for OVI emission as a function of red- 
shift. The surface brightness PDFs are similar at low fluxes, but 
diverge at the highest fluxes. In particular, the highest fluxes are 
found at the lowest redshift, which may be mostly due to the 
varying physical size of the pixels. 



6 WHAT TYPE OF GAS DOMINATES THE 
EMISSION? 

In this Section we investigate which gas contributes most 
of the UV emission. In particular, we study the dependence 
of the emission on density and temperature in Section [6.11 
We present the correlations between the observed fluxes and 
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Figure 7. Comparison of tlie intensities of the lines listed in Table [T] The contours show the number density of pixels and are logarith- 
mically spaced by 1.5 dex. The distributions have been calculated using five 20 comoving Mpc thick slices through the simulation 
box at ^ = 0.25, assuming an angular resolution of 15". Both axes show the line surface brightness Logjg5B in units of photon s~^ cm~^ 
sr~^. The dashed diagonal line in each panel indicates equal fiuxes. The CIII and CIV lines produce the highest fiuxes. Lines whose 
emissivities peak at similar temperatures (i.e. panels near the diagonal running from the bottom-left to the top-right of the plot) are 
strongly correlated. 



the physical properties of the gas in Section 16.21 and we 
finally calculate the emission from gas in different density 
and temperature ranges in Section [6.31 



6.1 Temperature-density plane 

In Paper I we demonstrated that aX z — 0.25 the bulk of the 
metals do not trace the bulk of the IGM mass. In fact, most 
gas mass resides in the low density IGM, while metals are 
more common in the m oderately dense, WHIM gas (see also 
|w lersma et al.ll2009bl ) and in the relatively dense gas with 
riH > 0.01 cm"'' and T < 10^ K that is cooling on to galax- 



ies. As stated in Section [2.31 we neglect the emission from 
the interstellar medium, assumed here as all the gas with 
riH > 0.1 cm~^, because the simulations lack the resolution 
and the physics to model its multi-phase composition. We 
found in Paper I that the X-ray emission traces neither the 
mass nor the metals in the IGM, but the densest and most 
metal enriched phases with temperatures close to the peak 
temperatures of the emissivity curves. 

In Fig. [To] we show the emission- weighted temperature- 
density distributions of the gas for our sample of UV lines, 
where the weight is given by the emission of a different line in 
each panel. The vertical line at wh =0.1 cm""' indicates the 
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Figure 10. The UV emission-weighted temperature-density distribution at z = 0.25. The distribution is weighted by the emission of a 
different ion, as indicated in the top left corner of each panel. The ions we consider here are CIII, CIV, Si IV, NV, OVI and NeVIII. 
The contours are logarithmically spaced by 1.5 dex and represent equal emission levels. The horizontal lines at constant temperature 
highlight the WHIM range, while the vertical line at nu = 0.1 cm"'' indicates the limit above which we impose an effective equation 
of state on to star-forming particles. For reference, the cosmic mean density corresponds to njj ~ 10~^ cm"''. CIV emission mostly 
traces dense warm gas in the haloes of galaxies. The bulk of the emission from O VI arises from equally dense, but hotter gas, while N V 
emission traces an intermediate phase of gas between those of CIV and OVI. NeVIII emission is produced by hot gas with T > 10^ K, 
whose distribution partially overlaps with that of OVI, but is only marginally traced by CIV lines. 



density threshold above which star forming particles obey an 
effective equation of state. The horizontal lines at T = 10^ 
K and 10^ K delimit the temperature range of the WHIM. 
We find that most of the emission traces moderately over- 
dense gas with nn > 10"''' cm~^ and that the bulk of the 
emission clusters around the emissivity peak temperature of 
each line. As such, NeVIII is a good tracer of hot gas with 
T ~ 10*^ K, while N V and OVI trace the cooler fraction of 
the WHIM with temperatures in the range 10^ K< T < 10® 
K. C IV traces gas slightly cooler than the WHIM, with 10'' 
K< T < lO^-'^ K. The bulk of the CIII and Si IV emission 
arises from gas with 10^ K< T < 10* K. The large scatter in 
some of the distributions shows that small fractions of the 
C IV, N V and O VI emission trace gas with temperatures 
that can differ by more than an order of magnitude from 
the peak temperatures of the emissivity curves, while the 
relatively narrow distributions in temperature of the C III 
and Si IV emissions indicates very inefficient emission out- 
side that narrow temperature range. 



6.2 Correlations between luminosity and gas 
properties 

The distributions in Fig. [10] show what gas produces most 
of the emission in each line. However, a large fraction of 
the energy is emitted at low fiuxes that are impossible to 



observe in the foreseeable future. It is therefore interesting 
to investigate the typical properties of the gas that produces 
fiuxes that might be within reach of future observations. 

We do this in Fig. 1111 which shows the median den- 
sity (left panel), temperature (middle panel) and metallic- 
ity (right panel) of the g function of the emitted fiux. 

The scatter in the distributions is large at low fiuxes, and 
decreases with increasing fiux. The flux on the x-axis is the 
fiux emitted by individual particles and the density, temper- 
ature and metallicity are local quantities which cannot be 
compared directly with the global gas properties that can 
be inferred observationally for the ICM and IGM. This is 
particularly relevant for the gas metallicity, because in the 
simulations the metal distribution is very inhomogeneous on 
small scales. 

As observed in Paper I, the median density and metal- 
licity of the gas increase steadily with flux, as would be ex- 
pected from the Zp^ dependence of the emissivity. On the 
other hand, the median gas temperature varies strongly for 
the lowest flux, then flattens out at the highest fluxes with 
F > 10"^° photon s"^ cm . This confirms the finding of 
Fig. [To] that most of the emission comes from gas with tem- 
peratures close to the peak temperature of the emissivity 
curve. Indeed, the highest fiuxes are produced by gas in a 
narrow temperature range centred on this value. As such, 
the detection of emission lines gives stringent constraints on 
the gas temperature and provides information on both the 
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Figure 11. Median values of the particle gas density (left panel), temperature (middle panel) and metallicity (right panel) that produce 
a given UV particle flux at 2 = 0.25 for the lines CIII, CIV, Si IV, N V, O VI and NeVIII. Note that fluxes are given in units of photon 
cm~^, but that previous figures used photon s~^ cm~^ sr~^. For reference, pmcan(2 = 0.25) corresponds to rtjj si 4 X 10~'^cm~'^ 
and our star formation threshold corresponds to logjQ p/pmean ~ 5.4. Emission from higher density gas (i.e. the ISM) was ignored. The 
median density and metallicity increase with luminosity, while the median temperature flattens out at the temperature for which the 
emissivity of each line peaks. 
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Figure 12. As Fig. O but showing maps of O VI emission from gas in four different temperature ranges. From left to right: T < 10^ K, 
IQS K < T < 10^ K, 10^ K < T < 10^ K and T > 10^ K. The highest O VI fluxes are produced by gas with 10^ K < T < 10^ K. 



density and metallicity of the gas. CIII and Si IV lines are 
the only exceptions to this rule, and show a slight increase 
in the median temperature at the highest fluxes, while never 
fully reaching a plateau. The detection of multiple lines may 
further help to break the degeneracy between density and 
metallicity. 

We have verified that the trends with density and tem- 
perature persist if metallicity is assumed to be constant. 

6.3 Density and temperature cuts 

In this Section we create O VI emission maps using only gas 
with temperatures and densities in a given interval. This 
allows us to probe the dependence of the emission on tem- 
perature and density and to visualise how the emitting gas 
is distributed in space. 

Fig. 1 121 shows O VI emission maps made using only gas 
particles that have temperatures in a specific range. We 



consider the four temperature intervals T < K, 10^ K 

< T < 10^ K, 10® K < T < 10'' K and r > 10^ K. Similarly, 
in Fig. [13] we consider four different density intervals, with 
LogioP/pmoan Varying by one dex in each bin. The surface 
brightness PDFs for the same cuts are shown in Fig. [14] for 
the temperature cuts and in Fig. [TS]for the density cuts. 

Figs. [T2lil5l demonstrate that the strongest OVI emis- 
sion comes from gas in dense (Logj^gp/pmcan > 2) and mod- 
erately hot (lO'^ K < T < 10® K) regions. Gas with 10® K 

< T < 10^ K contributes a significant fraction of the bright 
emission, while cooler gas with T < 10^ K contributes only 
a minor fraction of the total emitted energy. The O VI emis- 
sion of hot gas with with T > 10^ K is negligible. Compar- 
ing Fig. [12] with Fig. [13] shows that the warm gas emitting 
in O VI can be identified with the circum-galactic gas. Ob- 
serving the diffuse filamentary gas (lO'' K < T < 10® K; 
1 < p/pmean < 2) would require detection of emission with 
5*3 ~ 1 photon s^^ cm~^ sr~^. 
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Figure 14. As Fig. [5] but siiowing tiie PDF of ttie O VI surface 
briglitness from gas in different temperature intervals. The higliest 
OVI flux is produced by gas witli 10^ K < T < 10^ K, wtiile 
tlie emission is weakest in the highest temperature range. The 
detailed flux distributions vs. temperature vary line by line, and 
strongly depend on the peak of the emissivity curve of each line, 
with the strongest fluxes associated with the temperature range 
where the emissivity peaks. 

7 DEPENDENCE ON THE PHYSICAL 
MODELS 

The principal feature of the OWLS project is the availabil- 
ity of a large number of runs with varying physical prescrip- 
tions. In this Section we investigate how a number of phys- 
ical implementations affect the predictions for the WHIM 
emission. We present results for OVI A1032 A emission, 
but note that the same qualitative conclusions hold for the 
other lines considered in this study. We focus on models with 
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Log,o Sg (photon s"' cm"^ sr"') 

Figure 15. As Fig. |5] but showing the PDF of the OVI sur- 
face brightness from gas in different density intervals. The high- 
est fluxes are associated with the densest gas. In particular, the 
maximum surface brightness predicted for gas outside the virial 
radii of haloes (i.e. p/ pmcan J; 200) is at least a few orders of mag- 
nitude weaker than the maximum surface brightness from denser 
regions. 

varying feedback schemes and cooling recipes, which are the 
physical implementations that produce the largest variations 
for the diffuse emission. Physical prescriptions are changed 
one at a time, in order to isolate their impact on the results. 
The physical model s considered here ar e the following (we 
refer the reader to ISchave et al] (|2010l ) and to the works 
cited below for a more detailed description of the specific 
physical implementations) : 

(i) REF: Our reference model, as described in Section [2TT] 
All results presented in the preceding sections were obtained 
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Figure 16. As Fig.[3l but showing O VI emission maps for different simulations, as indicated in the top left eorner of each panel. 



from this si mulation. This model uses the st ar formation pre- 
scription of ISchave fc Dalla Vecchial (|2008|). th e sup ernova- 
driven wind model of Dalla Vecchia fc "Schavel jiooa ) with 
wind mass loading 77 = 2 (as defined in Section 12.11) and 



wind initial velo city « 
cooling rates of 
lution model of 



Wiersma et al 



600 k m s 

2OO9J 



Wiersma et al 



, the metal-dependent 
the chemical evo- 



( 2009bl ) and the WMAP3 



cosmology. AGN feedback is not included. 

(ii) NOSN: As REF, but without supernova feedback. 
The metals produced by stars are only transported by gas 
mixing and no energy is transferred to the IGM and ISM 
when SNe explode. This model is clearly unrealistic, it 
for example predict s far more star formation than obseved 
l|Schave et al.ll2010l '). We include it, however, because it can 



shed light on the physics that determines the metal-line 
emission. 

(iii) NOZCOOL: As REF, but with radiative cooling rates 
calculated assuming primordial abundances. The calculation 
of the metal line emission in this model is not self-consistent, 
since the contribution of metals to the cooling rate is not 
included. The neglect of the contributions of metal lines im- 
plies slower gas cooling than in the REF model. We in- 
clude this model to investigate the role of metal-line cooling 
and because previous work often a ssumed primordial abun- 
dances (e.g. iFurlanetto et al.ll20o3 . see Section[8]). 

(iv) AGN: As RE F, but with the additio n of AGN feed- 
back as described in Booth fc Schavel (|2009l) . We note that 
[McCarthy et al.l (|20ia )" have shown that this model (but not 
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Figure 17. Same as Fig. 1161 but for constant gas metallicity Z = O.IZq. 

600 km and 770 = 2. This implies that at pcrit the gas 
is kicked with the same velocity and mass-loading as in the 
REF model. This particular scaling of wind velocity with 
density was chosen so that the wind velocity scales with the 
local sound speed (recall that we impose the equation of 
state P oc p*/"^ onto the ISM) and the wind mass-loading 
is scaled such that the total wind energy is independent of 
the gas density. As a consequence, the wind energy per unit 
stellar mass is equal to that in the REF model. 

(vii) WVCIRC: As REF, but with "momentum- 
driven" galactic winds f ollowing the prescription of 
lOppenheimer fc Pavel (|2008t ) , with the difference that winds 
are "local" to the star formation event and are fully 
hydrodynamically coupled as in the REF model (see 



model REF) reproduces the observed density, temperature, 
entropy and metallicity profiles of groups of galaxies, as well 
as their star fractions. Model AGN may therefore, for our 
purposes, be the most realistic model of the ones considered 
here. 

(v) WML4 : As REF, but with a wind mass loading twice 
as high, that is 77 = 4. The total energy injected in winds is 
also twice that in REF. 

(vi) WDENS: As REF, but the wind properties scale with 
the local gas density as: 

= V^o (p/pcrit)^''^ (3) 

r) = 770 (p/PcHt)-'/^ (4) 
where pcrit is the density threshold for star formation, v^o = 
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Figure 18. As Fig.[5l but showing the surface brightness PDFs of O VI emission for simulations with varying physical prescriptions. The 
simulations are the same as in Fig. 1161 The results are surprisingly robust to changes in the physical model. However, neglecting metal-line 
cooling results in an overestimate of the flux at the bright end. Ignoring supernova-driven winds has little effect for the brightest pixels, 
but otherwise shifts the PDF to lower fluxes. AGN feedback significantly reduces the maximum flux. 




Figure 19. As Fig. 1181 but for simulations with varying physical proscriptions, under the assumption of a constant gas metallicity 
of Z = 0.1^0. The simulations shown are the same ones as in Fig. 1161 The assumption of constant metallicity boosts the flux from 
the low-density, cool regions of the universe, where the largest contribution to the line emissivity comes from photo-ionisation by the 
UV background (compare the inset to that of Fig. list , but it reduces the maximum fluxes. The differences between the models are 
qualitatively the same as when the predicted metal distributions are used (compare with Fig. I18II . 
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iDalla Vecchia fc Schavell2008l for a discussion of the impor- 
tance of these effects). The wind initial velocity and the 
wind mass loading factor are defined as a function of the 
galaxy velocity dispersion a cLS — 5(T and iq = i>wo/o- 
respectively, with «wO = 150 km s~^. The velocity disper- 
sion a = v^fcirc is estimated using an on-the-fiy friends-of- 
friends halo finder. This model was motivated by the idea 
that galactic winds may be driven by radiation pressure on 
dust grains rather than SNe. The energy injected into the 
wind becomes much higher than for REF for halo masses 
exceeding lO" - 10^^ Mq. 

(viii) MILL: As REF, bu t using the WMAP year 1 cos- 
mology jSpergel et al.ll20o3 ) as in the Millennium simulation 
l|Springel et al.ll2005l ) and a wind mass loading factor = 4. 
Results from this model must be compared to those of the 
WML4 model to isolate the effect of cosmology. 

(ix) DBLIMF: As REF, but using a top-heavy IMF at 
high gas pressures. This model switches from a Chabrier 
IMF to a power-law IMF dN/dM oc (as compared to 
oc M~^-^ for the high-mass tail of the Chabrier IMF) at 
the critical pressure P/k — 2.0 x lO'' cm~'^ K, which was 
chosen because for this value ~ lO"'^ of the stellar mass 
forms at higher pressures. The extra SN energy is used to 
increase the initial wind velocity to 1618 km s~^for star 
particles formed with a top-heavy IMF, while keeping the 
mass loading factor fixed at = 2 (as for the Chabrier IMF, 
these parameter values correspond to 40% of the available 
SN energy). This model is called DBLIMFCONTSFV1618 
in the OWLS project (which i ncludes more model s with top- 
heavy IMFs in starbursts, see lSchave et al.ll201Cll l. 

Fig. ll6l shows maps of the O VI emission for all runs in a 
region of 14 h^^ Mpc on a side, or equivalently 1.12 degrees 
on the sky. Fig. [T7] shows the emission in the same simula- 
tions, but assuming a constant gas metallicity oi Z = Q.IZq. 
The corresponding surface brightness PDFs of the OVI 
emission are presented in Fig. 1181 while those for constant 
metal abundances are shown in Fig. 1191 

Despite the substantial differences in the physical mod- 
els, we find that most models produce relatively small differ- 
ences in the results. In particular, the specific implementa- 
tion of stellar feedback ( WVCIRC, WDENS) barely affects 
either the surface brightness PDF (Fig. 1181) . Even the small 
differences between feedback models found in Paper I for 
the X-ray emission are mostly washed out when considering 
lower ionisation transitions in the UV band. The only excep- 
tions that show noticeable differences with the REF model 
are NOSN, NOZCOOL and AGN and we now discuss each 
of these in turn. 

The spatial distribution of the emission for the NOSN 
model is more clumpy and concentrated around galaxy 
haloes than for all other models. The fiux PDF of the NOSN 
model differs primarily at low fiuxes (< 1 photon s^^ cm^^ 
sr~^). This is a consequence of the fact that without SN 
winds the transport of metals to the IGM is very ineffi- 
cient. This naturally produces very low IGM metallicities, 
and therefore low metal line emission away from galaxies. Al- 
though the surface brightness distribution predicted by the 
NOSN model differs greatly from that in the runs that in- 
clude SN feedback (Fig. [16]), when the same maps are made 
assuming constant gas metallicity (Fig. I17p most of the dif- 
ferences disappear. This indicates that the primary way in 



which SN feedback affects the UV emission is through alter- 
ing the metal distribution and not necessarily by changing 
the thermodynamic properties of the gas emitting in the UV 
band. 

The NOZCOOL run predicts a larger fraction of pixels 
with high fluxes than all other models. This indicates that 
the NOZCOOL model predicts a larger fraction of gas within 
the temperature range that produces O VI emission. 

The surface brightness PDF for the AGN model is re- 
markably similar to those of other models for all but the 
highest fluxes (> 100 photon s~^ cm~'^ sr~^), where it shows 
a steep cut-off. This indicates that the flux from high den- 
sity regions is lower than in the other models and is likely 
to be a consequence of two physical processes. Firstly, AGN 
feedback suppresses star formation in haloes that host black 
holes, which in time would substantially reduce the gas 
metallicity. Secondly, a s demo nstrated by iMcCarthv et al.l 
|2010| ) and iDuffv etaP (j2010l ). the inclusion of AGN feed- 
back in the simulations significantly decreases the gas frac- 
tions in the centres of all massive haloes. Together, decreases 
in both the density and metallicity of the gas lead to greatly 
suppressed fiuxes in the highest density regions. 

Variations in the cosmology, seen by comparing the 
MILL and WML44: models, do not produce significant vari- 
ations in the PDF, but do affect the spatial distribution of 
the emission, as can be seen in Fig. 1161 This is mostly a con- 
sequence of the fact that structure formation has progressed 
further in the MILL run due to the higher value of erg used 
in this simulation. 

If we assume a constant gas metallicity, as in Figs. [17] 
and 1191 then the results barely vary with varying physical 
prescriptions, with the exception of the NOZCOOL model, 
whose results reflect those for the case with the simulation 
metallicities. 

In summary, except for metal-line cooling and AGN 
feedback, the physical prescriptions we have tested here only 
marginally affect the properties of the gas reservoir that pro- 
duces UV emission, which is therefore insensitive to varia- 
tions in the SN feedback scheme and stellar IMF (as in the 
DBLIMF model). The results presented in this work are 
therefore robust against uncertainties in the subgrid mod- 
els discussed in this section, although AGN feedback may 
significantly suppress the peak fiuxes. 



8 COMPARISON WITH PREVIOUS WORK 

Our predictions for C IV and O VI emission can be com- 
pared directly to those of Furlanetto et al. (2004, F04 
hereafter, Figs. 7-10 ), who use the "G5" simulations of 
ISpringel fc Hernguistl l|2003l '). In particular. Fig. 10 of F04 
uses the same box size and redshift as our Figs. [5] and 1181 
F04 use a slice thickness of Az — 0.01, whereas Az « 0.0074 
for us. Since the bright-end of the surface brightness PDF 
is proportional to the slice thickness (see appendix I A3|) . we 
have reduced their PDFs by about 0.13 dex for a fair com- 
parison. They use a pixel size of 20 proper kpc//i, which 
is about half the size of our pixels. However, Fig. [8] shows 
that this difference has little effect on the results. The sim- 
ulation of F04 uses about four times less particles, but our 
Fig. [XT] indicates that this is sufficient. They use a WMAPl 
cosmology instead of WMAP3, but Fig. 1181 shows that this 
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also does not make a significant difTerence. The simulation 
used by F04 used different prescriptions for chemodynamics 
(they did not track individual elements and they did not 
follow the delayed release of heavy elements by stars), star 
formation and SN feedback and, like our reference model, 
did not include AGN feedback. The most important differ- 
ence is, however, likely the fact that F04 did not include 
metal-line cooling. 

Comparing F04's prediction for O VI to that for our 
REF model, we find that bright regions are somewhat more 
common in their simulation. However, since they did not 
include metal-line cooling, it is more appropriate to compare 
with our model NOZCOOL. Reassuringly, in that case the 
bright ends of the PDF agree nearly exactly. The same, is 
not true, however, for C IV. While F04 predict the maximum 
fluxes to be weaker for C IV than for O VI, we find the 
opposite. Since our predictions agree for O VI (if we ignore 
metal-line cooling, which is not self-consistent), this means 
we predict substantially more C IV emission. This could be 
because F04 did not include metal enrichment by SN la and 
AGB stars. While this is not that important for oxygen, 
which is released mostly by massive stars, intermediate mass 
stars are important for carbon. 



9 CAN THE WHIM BE DETECTED IN THE 
UV? 



In Section 16.31 we showed that UV emission lines such 
as CIV, NV, OVI and NeVIII are good tracers of the 
coolest fraction of the WHIM and tend to be strongest 
in the haloes of galaxies and outside the central cores 
of groups. In this Section, we discuss their observability 
by current and future instruments, in particular by COS 
iFroning fc GreenI 120081'), FIREBAL L jTuttle et al.l I2OO8I ') 
and ATLAST ( |Postman et al.ll2008l 'l. The requirements for 
UV telescopes to detect metal-line emission lines from the 
WHIM are the same as those for X-ray instruments: a large 
field of view (FoV) plus high spatial and angular resolution. 

COS, which was successfully mounted on the Hubble 
Space Telescope, has been designed to detect WHIM ab- 
sorption features, rather than diffuse emission. Its extremely 
small FoV is optimal to detect light from faint UV point 
sources, but its spectral resolution declines steeply for ex- 
tended sources. The flux limit of COS, of order 1000 pho- 
ton s~^ cm~'^ sr~^, is too high for allowing the detection of 
WHIM emission in filaments or outside groups at z — 0.25, 
which is of the order of 0.01 — 1 photon s~^ cm~^ sr~^. How- 
ever, it might just be sufficient to detect metal-line emission 

from the densest gas in the haloes of galaxies. 

The FIREBALL balloon experiment jTuttle et al.1 

I2OO8I ) has been successfully launched in June 2009, with 
the specific aim of detecting diffuse emission from C IV at 
z ~ 0.3 and OVI at z ~ 1. FIREBALL is a pathfinder inte- 
gral field unit with a relatively large field of view of 2.5' x 2.5' 
and an angular resolution of 10" . The spectrograph is fed by 
about 400 fibres sensitive in the narrow UV band 1950A- 
2200 A, and can reach a fiux limit of ~ 2000 photon s cm 
sr~^. Such a fiux limit is substantially lower than the max- 
imum value of ~ 10^ photon s~^ cm~^ sr~^ we predict for 
CIV emission at z ^ 0.25 (see Fig. [5]). For OVI we predict 
a maximum surface brightness at z ~ 1 of ~ 3000 photon 



s~^ cm~^ sr~^ (see Fig. [9]), which is close to the detection 
threshold. However, AGN feedback may decrease these max- 
imum fiuxes by an order of magnitude (Fig. I18|) . As such, 
we predict that FIREBALL should be capable of detecting 
a significant number of photons associated with diffuse C IV 
emission from high density regions. Detection of O VI emis- 
sion from high density regions might also be feasible, but 
this represents a more challenging prospect at the very lim- 
its of the instrument capabilities. C III is predicted to be 
even stronger than C IV, but this line does not trace the 
WHIM and, being a singlet, it would be dif ficult to separate 
from contaminating HI Ly series lines (e.g. lFurlanetto et al.l 
l2005h . We are looking forward to seeing the first results. A 
positive detection of diffuse emission would certainly be a 
powerful argument in favour of future larger UV telescopes. 

Ideally, next generation UV missions targeted to detect 
diffuse emission from the WHIM should aim to achieve flux 
limits of the order of 1-10 photon s~^ cm~^ sr~^ within 
a FoV with size of at least an arcminute. Such a low flux 
limit is necessary to detect emission lines from truly diffuse 
gas outside the haloes of galaxi es. The proposed spe cifica- 
tions for the ATLAST mission l|Postman et al.1 1200^ 1 may 
be ideally suited for the detection of WHIM emission from 
fllaments. 

Finally, we note that we have not considered emission 
from hydrogen Lyman series lines, because they are thought 
to trace cooler gas than the WHIM. However, the intensity 
of LyQ emission may well be higher than that of the metal 
lines and Lyman series lines are therefore a potential con- 
taminant fo r other rest-frame UV lines (e.g. lFurlanetto et al.l 
l2003l . l2005l ). Although doublets can be unambiguously iden- 
tified with sufficient spectral resolution, singlets like C III 
(A977 A) would be difficult to distinguish from contaminat- 
ing hydrogen lines unless other lines corresponding to the 
same redshift are detected. 



10 CONCLUSIONS 

UV metal lines constitute one of the main cooling channels 
of warm gas. These UV lines provide an attractive route 
towards detecting the cooler fraction of the missing baryons 
in the low redshift Universe, which simulations predict to 
reside in a warm- hot diffuse phase with 10^ K < T < lO'' K 
(the WHIM). 

In this work, we used a subset of cosmological simula- 
tions fro m the Overwhelmin gly Large Simulations (OWLS) 
project dSchave et al.l I2OIOI ) to study the nature and de- 
tectability of rest-frame UV metal-line emission from dif- 
fuse gas (nu < lO-^cm-3). The OWLS runs are among 
the largest dissipative simulations ever performed, and use 
new prescriptions for a range of baryonic processes, such 
as radiative cooling, star formation, and feedback from SNe 
and AGN. Of particular relevance is that the OWLS runs 
are fully chemodynamical and track the release of 11 differ- 
ent elements by massive stars, SNe of types II and la, and 
asymptotic giant branch stars. In addition, radiative cool- 
ing is implemented element-by-element in the presence of an 
evolving UV/X-ray background. 

Our conclusions can be summarised as follows: 

(i) In the UV band the strongest metal lines are C III 
(977 A, T ~ 10^ K), CIV (1548,1551 A, T ~ 10^ K) and 
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O VI (1032,1038 A, T ~ 10^'^ K). C III is the strongest line 
and mostly traces circum- galactic gas. Si IV (1393,1403 A) 
traces a fraction of the high density gas traced by C III emis- 
sion, within a smaller range around the peak temperature 
T ~ 10^ K. The maximum C IV emission is almost an or- 
der of magnitude weaker than the peak C III emission, and 
traces gas that has a similar spatial distribution as the gas 
that produces CIII emission, but is slightly warmer. OVI 
emission arises from warmer and more diffuse gas than C IV 
emission. As a consequence, O VI emission is weaker than 
C IV close to galaxies, but stronger outside haloes. The max- 
imum surface brightness for N V (1239,1243 A, T ~ lO^ ''^ K) 
and Ne VIII (770,780 A, T ~ 10*^ K) is weaker than that of 
O VI by up to an order of magnitude. The properties of N V 
emission are intermediate between those of C IV and O VI. 
The NeVIII line is qualitatively similar to OVI and traces 
warm-hot gas with T ~ l(f K. 

(ii) Emission lines provide strong constraints on the tem- 
perature of the emitting gas, because the bright emission 
traces gas with temperatures close to the peak of the emis- 
sivity curve. Photo-ionisation by the UV background is 
unimportant for the gas that dominates the emission from 
the intergalactic medium. Since this temperature increases 
with atomic number, C IV lines trace cooler gas than O VI 
emission, which in turn comes from cooler gas than Ne VIII 
emission. 

(iii) UV emission lines are likely to be detectable only in 
relatively high density regions. The detection of gas with 
density p < lO^pmcan requires Sb, ^ 10^ photon s~^ cm~^ 
sr~^, while gas with p < lO^pmcan requires Sb, ^ 10 photon 
s~^ cm~'^ sr~^. The gas that would be detected at these 
surface brightness levels tends to be relatively metal-rich 
{Z > 10~^ Zq ). The detection of low density structures such 
as filaments (p ~ lOpmoan) is thus extremely challenging 
and might be best achieved through absorptio n techniques. 
Howe ver, the FIREBALL balloon experiment (|Tuttle et al.l 
|2008| ) has a good chance to detect diffuse C IV emission at 
z ~ 0.3. We will present comparisons with absorption line 
data elsewhere 

(iv) By comparing results from different OWLS models, 
we have investigated the dependence of the results on a num- 
ber of physical mechanisms. The variations in the predicted 
fluxes for different simulations are within a factor of ten. 
Varying the intensity of feedback or the stellar IMF affects 
the simulated fluxes by at most a factor of a few. Assum- 
ing primordial abundances for the calculation of the cooling 
rates (which is not self-consistent, but often done in previ- 
ous work) boosts the flux in the brightest pixels by about an 
order of magnitude. Simulations that include AGN feedback 
predict maximum fluxes about an order of magnitude lower 
than the REF model in massive haloes, mostly because of 
the lower gas metallicity that results from the suppression of 
star formation. However, AGN feedback has little effect on 
the surface brightness PDF for values more than an order 
of magnitude below the maximum value (i.e. for Sb ^ 10^ 
photon s~^ cm~'^ sr~^). 

(v) The quadratic dependence of the emissivity on density 
implies that emission is a highly biased tracer of the IGM 
mass and metals. While emission lines are thus not an ideal 
tool to close the baryon budget, they are excellent tracers of 
the gas flowing in and out of galaxies. Rest-frame UV (and 
soft X-ray) lines are sensitive to the gas temperatures that 



are relevant for galactic winds, cooling flows and cold-mode 
accretion. Thus, they have the potential to open up a new 
window onto some of the most poorly understood aspects of 
the formation and evolution of galaxies. 

By comparing the results presented in this work with 
those of Paper I, we find that, on average, UV emission lines 
have a surface brightness two to three orders of magnitude 
higher than soft X-ray lines. This can be seen, for example, 
when comparing O VI and O VIII emission. The maximum 
surface brightness predicted for O VI, assuming an angular 
resolution of 15", is about 10^ photon s~^ cm~^ sr~^, while 
for O VIII it is ~ 10^ photon s"^ cm ^ sr ^, or two orders 
of magnitude lower. It has to be noted, however, that the 
relative strength of the emission lines varies strongly with 
the gas temperature. As a consequence. X-ray and UV lines 
are well suited to probe gas with different properties that 
may not physically occupy the same regions of space. 

Effectively, UV and soft X-ray metal lines together 
have the potential to yield a huge amount of information 
about the WHIM and the circum-galactic medium. While 
X-ray lines provide clues about the hottest component with 
T > 10^ K, UV lines trace the cooler fraction of the gas 
with T < lO" K including WHIM, galactic winds, cool- 
ing flows and cold mode accretion. Ultimately, the multi- 
wavelength approach may constrain the density, tempera- 
ture, metallicity and ionisation state of gas, highlighting the 
fully 3-dimensional structure of the gas distribution. This 
would shed new light on a wealth of processes, from the 
feedback mechanisms responsible for the metal enrichment 
of the IGM, to the chemical evolution of galaxies, and to the 
cosmic flows that provide fresh fuel for star formation and 
black hole accretion. 
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APPENDIX A: CONVERGENCE TESTS 

In this Appendix we investigate how a number of numerical 
issues affect our results. In particular, in Appendix lAll we 
vary the mass resolution of the simulations, while in Ap- 
pendix IA2I we keep the mass resolution constant and vary 
the size of the simulated volume. Finally, in Appendix I A3I 
we show how varying the slice thickness affects the shape of 
the surface brightness PDF. All tests use the REF physical 
model and results are shown only for the O VI 1032 A line. 

Al Mass resolution 

The effect of varying the numerical resolution of the simu- 
lations is shown in Fig. lAll which gives surface brightness 
PDFs for maps created for three different realisations of the 
same simulation with particle numbers A'' = 2x 128'^, 2x256'^ 
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Figure Al. As Fig. [S] but showing the O VI surface brightness 
PDF as a function of the resolution of the simulation. The particle 
mass (gravitational softening) increases by a factor of eight (two) 
when the number of particles decreases by a factor of 2^ . All runs 
assume the same physical model, angular resolution and box size. 
The predictions of the A'^ = 256"^ and 512'^ runs are very close 
for fluxes greater than 10~^ photon s~^ cm~^ sr~^, but diverge 
again for the small number of pixels with fluxes greater than 
10'^ photon s~^ cm~^ sr~^. 

and the default number of 2 x 512^ . The dark matter particle 
masses are Mdm = 2.6 x 10^° h'^ Mq, 3.2 x 10^ h'^ Mq 
and 4.1 X 10* M!q, respectively. All runs have a box size 
of 100 Mpc and we assume an angular resolution of 15". 

The results for the two high resolution runs with A'' = 
2 X 256^ and 2 x 512"^ particles converge well for intermedi- 
ate fluxes, but diverge somewhat at the highest and lowest 
fluxes, with differences within a factor of a few. On the con- 
trary, the lowest resolution run with N = 2 x 128^ does not 
converge anywhere and the overall shape of the PDF is dif- 
ferent from those of the other runs, with a lot more power at 
the lowest fluxes. This is due to the fact that this simulation 
is not able to properly describe the density and temperature 
distributions of the gas both on large and on small scales. 
The star formation history is severely underestimated in the 
TV = 2 X 128^ simulation and the metals produced by stars 
are not efficiently transported into the diffuse IGM. This im- 
plies a far lower gas metallicity than in the other runs and 
explains the very low fluxes observed outside haloes. 

A2 Box size 

The effect of varying the box size of the simulation on the 
O VI surface brightness PDF is shown in Fig. IA2I The runs 
assume different initial conditions and allow us to test how 
the large scale structure affects the flux statistics. The mass 
resolution in the simulations is kept constant, with a dark 
matter particle mass Mdm = 4.1 x 10® /i"^ Mq. 

We flnd that the surface brightness PDFs converge, but 



Figure A2. As Fig. [S] but showing the O VI surface brightness 
PDF for simulations with box sizes of 100, 50 and 25 h''^ Mpc. All 
runs assume the same physical model, particle mass and angular 
resolution. Except for very low fluxes (< 10~^ photon s~^ cm~-^ 
sr~^), the results are insensitive to the size of the simulation box. 
The highest flux values that arc sampled increase with box size. 

that the highest fluxes are only sampled in the largest box, 
in which more extreme objects can form. The surface bright- 
ness PDF for the 25 ft"^ Mpc run diverges from the other 
two at the lowest fluxes. 



A3 Slice thickness 

In this Section we investigate the effect of varying the thick- 
ness of the slices we cut through the simulation box on the 
shape of the surface brightness PDF. This allows us to test if 
the strongest emission, coming from haloes of different sizes, 
is entirely contained within a slice. Results are shown in Fig. 
IA3l for four different slice thicknesses (40, 20, 10 and 5 
Mpc). To facilitate the comparison, we have normalised the 
surface brightness PDF by the slice thickness Aa;. 

Within the potentially detectable regime the results are 
insensitive to the slice thickness. This occurs because the 
slices considered here are thick enough to fully contain the 
collapsed objects that give rise to the bulk of the UV emis- 
sion, but at the same time they are thin enough that they 
avoid the superposition of multiple structures, which would 
strongly increase the flux in those pixels where more than 
one structure could be found along the line of sight. At lower 
fluxes, the shape of the PDF depends on the flux in low den- 
sity regions and on the number of low density structures such 
as fllaments that can be found in each slice, which depends 
on the slice thickness. As a consequence, the width of the 
PDF at intermediate fluxes scales with the inverse of the 
slice thickness and the number of pixels at the lowest fluxes 
increases for decreasing slice thickness. 

Throughout this work, we have assumed a slice thick- 
ness of 20 Mpc. Fig. IA3I demonstrate that this choice 
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Figure A3. The OVI surface brightness PDFs, normahsed by 
the slice thickness Ax, for maps with varying slice thickness: 40, 
20, 10 and 5 Mpc. The number of slices used to calculate 

each PDF is 2, 5, 10 and 20, respectively The test is for the REF 
model with box size 100 Mpc at z = 0.25. All maps assume 
the same angular resolution of 15". In the high flux regime the 
PDF is proportional to the slice thickness. 



is well justified and allows us to estimate the flux from high 
density regions accurately. For fluxes higher than about 1 
photon s~^ cm~^ sr~^ the PDF is proportional to the slice 
thickness. 

This paper has been typeset from a TJ^X/ I^TJ^X file prepared 
by the author. 
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